The approach to a superconductor-to-Bose-insulator transition in disordered films 
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Through a detailed study of scaling near the magnetic field-tuned superconductor-to-insulator 
transition in strongly disordered films, we find that results for a variety of materials can be collapsed 
onto a single phase diagram. The data display two clear branches, one with weak disorder and an 
intervening metallic phase, the other with strong disorder. Along the strongly disordered branch, the 
resistance at the critical point approaches Rq = h/4e 2 and the scaling of the resistance is consistent 
with quantum percolation, and therefore with the predictions of the dirty boson model. 
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At finite magnetic field, a disordered superconducting 
system, thinned down to a two-dimensional film, may ex- 
hibit a true superconducting state only at zero temper- 
ature. Upon increasing either the magnetic field or the 
disorder superconductivity will be lost Two main 
scenarios have been proposed for this transition. In the 
"dirty boson" scenario Cooper-pairs exhibit enough in- 
tegrity for the transition to be dominated by boson local- 
ization |2|. Alternatively, pairs are broken at the transi- 
tion and the system above the transition is dominated by 
fermion physics . In either scenario it is expected that 
above the transition a large enough system will exhibit 
insulating behavior as the temperature approaches zero. 
Focusing on the magnetic field-tuned transition, initial 
experiments seemed to confirm the boson-dominated sce- 
nario Q. However, scaling analysis around the appar- 
ent superconductor-insulator transition (SIT) [!, @, 0] 
yielded critical exponents that resemble those of clas- 
sical rather than quantum percolation, and the critical 
resistance at the transition was not the quantum of re- 
sistance for Cooper pairs, Rq = h/Ae 2 « 6.45 kfi/Q 
|8|. In addition, some experiments have challenged the 
general occurrence of a SIT, demonstrating an appar- 
ent transition that is only a crossover to a new metallic 
state at low temperatures. In all these cases the insu- 
lating phase above the transition was rather weak, sug- 
gesting the possibility that nearby fermions couple dissi- 
patively to the system, hence preventing the observation 
of a pure boson-dominated SIT [9(. More recent exper- 
iments on amorphous indium-oxide (InOx) |T3| and 
TiN fims [ll[ indicated that the strength of the insu- 
lating behavior above the apparent magnetic field-tuned 
SIT can be enhanced by using more disordered films 
However, no systematic study of the critical behavior was 
undertaken to shed light on the quest to achieve a true 
boson-dominated SIT. 

In this paper we study the SIT critical behavior of 
a wide range of InOx films at temperatures well below 
their mean-field superconducting transition (denoted as 
T o). While different batches of samples may have differ- 
ent properties for similar normal-state Ru 112} , we found 
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FIG. 1: Phase diagram of all samples including InOx and MoG e 
using reduced axes. Open circles (Ref. |2tl ) and squares (Ref. |1{J| ) 
are for MoGe films. Full circle is for thin Ta films (Ref. |21|). 
Diamonds are for InOx (Refs. [T7l . liaip and full triangles and squares 
are new InOx data for this work, from two batches. The error bars 
give a measure of the difficulty in determining H c 2(0). The solid 
line is a guide to the eye. The inset shows the low temperature part 
of the generic phase diagram as proposed in Ref. 0|. X denotes a 
possible critical point between a transition to a Bosc insulator and 
a transition to a metal (see text). The bottom shaded area denotes 
the region for which the conductance is smaller than h/e 2 . Samples 
"a" and "b" are discussed further in Fig. [2] Different regimes of 
the diagram as related to the data are explained in the text. 



that all samples can be collapsed onto a single phase di- 
agram if we plot the critical resistance at the SIT nor- 
malized by Rq versus the critical field at the transition 
normalized by H C 2(0) (i.e. the mean-field upper critical 
field at T = 0). Moreover, on the same phase diagram 
other samples from previously reported measurements on 
amorphous MoGe can be placed as well. Based on our 
observations we conclude that a true "dirty-boson" SIT 
can be achieved in the limit of strong disorder. In that 
limit the critical magnetic field approaches zero while the 
critical resistance approaches the quantum of resistance 
for pairs. Scaling near the critical point, as viewed from 
measurements at finite temperatures, is consistent with 
a quantum percolation critical point. In the opposite 
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limit of weak disorder, such as our earlier measurements 
on MoGe @, 0], the critical field approaches H C 2(0) and 
the critical resistance decreases with respect to Rq until 
it is indistinguishable from the normal-state resistance 
of the sample. In that limit scaling produces classical 
percolation critical behavior that is disrupted by the ap- 
pearance of an intervening metallic phase. Figure [1] is 
a summary of our result. It provides a new perspec- 
tive on the SIT quantum critical behavior, suggesting 
that while a true Bose-like SIT can be achieved in the 
limit of strong disorder, the other limit represented by 
the asymptote H C /H C 2(0) — 1 contains completely new 
physics for which a new quantum metallic state could be 
the featured highlight [13| . 

InOx films reported in this paper were fabricated by 
electron beam evaporation of a 99.999% pure indium ox- 
ide piece onto cleaned silicon nitride substrates [14] . The 
disorder was tuned during fabrication by controlling the 
sample thickness and the oxygen pressure [f| . The films 
were subsequently patterned into four-point resistivity 
bars by photolithography and argon-ion etching. Care 
was taken not to heat the samples above ~ 50°C in or- 
der to prevent any recrystallization [l5| . Resistance mea- 
surements were conducted in a dilution refrigerator using 
standard four-point lockin techniques. Over 15 samples 
were used for the present study. The mean-field transi- 
tion temperature T c q, an energy scale for the formation of 
Cooper pairs in these strongly fluctuating systems, was 
determined by fitting the fluctuation conductivity to the 
Aslamasov-Larkin model for two dimensions [l(| . Within 
each batch of samples disorder increases as T c o decreases 
and the normal resistance increases. The upper critical 
field was determined from the low field initial slope of 
T C (H) as explained in detail in 

Upon increasing the magnetic field, the behavior of the 
films changes from superconducting-like with resistance 
decreasing with decreasing temperature, to insulating- 
like with resistance increasing with decreasing temper- 
ature. The separating curve marks the SIT, defining a 
critical resistance R c and magnetic field H c . Isotherms 
of resistance versus magnetic field around this point all 
cross at H c with a critical resistance R c marking a zero- 
temperature quantum SIT phase transition. Assuming 
a magnetic field-tuned transition, the correlation length 
diverges at the transition as £,(H) oc \H — H c \~ v . The 
dynamics of the system are characterized by the relation 
between the correlation length and a vanishing frequency 
f2 through fi cx . At a finite temperature, Ml is cut-off 
by ksT, defining a length Lt for the cut-off of quantum 
fluctuations: ksT ~ L^ z - This leads to a scaling relation 
for the resistance which is commonly written as Q 
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Figure [2] shows the resistance isotherms in a perpen- 
dicular magnetic field for two InOx samples, one with 
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FIG. 2: Resistance isotherms for the two samples marked in 
Fig. [T] The isotherm temperatures are listed in the legend. Note 
the different scales of the vertical axis. Arrow points to the critical 
fields of H c = 12.6 T for sample "a", and H c = 0.83 T for sample 
"b" . Line A-B denotes the field at which the activation energy (as 
determined from the resistance vs. temperature at a fixed field) is 
maximum. 



weak disorder (Fig. [5^) and one with very strong disor- 
der (Fig. [2)3) . Similar types of samples were introduced 
earlier [5j. For the sample in Fig. [2^ the critical resistance 
is ~4.1 kil, very close to the normal state resistance be- 
fore superconductivity sets in, and the magnetoresistance 
peak is weak. Furthermore, the critical field is ~12.6 T, 
very close to H C 2(0) w 14 T calculated for this sample. 
Thus, this sample is very similar to the MoGe samples 
previously reported 0,0). Scaling of this sample, as well 
as all previously reported MoGe samples, typically fails 
at the lowest temperatures. However, in a finite temper- 
ature range scaling yields exponents similar to classical 
percolation (in particular v w 1.35). 

Starting from the crossing point and increasing the 
magnetic field towards the magnetoresistance peak we 
note in Fig. [5] that the resistance increases with decreas- 
ing temperature, at an increasing rate. At some field 
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this rate is maximum (denoted by lines (A-B)) and we 
can calculate an activation energy through this magnetic 
field cut (= Hp, also see [Hj). To evaluate the trend of the 
samples we also need a measure of the disorder. For each 
sample we take the maximum normal-state resistance at 
zero magnetic field, before superconductivity sets in, as 
such measure, denoted R™ ax . In Fig. [3] we plot H c and 
H p versus R™ ax , and we clearly identify three regimes: 
below the crossing point the phase is superconducting 
with resistance decreasing with decreasing temperature; 
at very high fields, above the field of maximum activation 
energy, the system is expected to be dominated by Fermi 
physics; in between the two sets of points the system be- 
haves as a Bose insulator The two solid lines, though 
drawn solely as guides to the eye, are not parallel but get 
farther apart as the disorder increases, indicating that the 
Bose-insulating phase is becoming more well-developed 
and less influenced by background fermions [f|. This is 
in contrast to figure 5 in Ref. [13] which uses a differ- 
ent measure for the sample disorder and is drawn over a 
smaller range of critical fields. 
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FIG. 3: Phase diagram for the superconducting and insulating 
phases. The various closed and open symbols denote the H c cross- 
ing and the H p activation peak, respectively, accounting for the 
two batches of samples. The lines are guides to the eye. 



InOx films appear to complement MoGe films in terms 
of the strength of the disorder, the two systems together 
spanning a wide range of the disorder axis. We now con- 
sider InOx samples for which the critical field H c is very 
low and the regime of Bose-insulating behavior is very 
large, and we analyze the critical behavior near the SIT. 
We have attempted to connect the results with materials 
of lower disorder such as MoGe, 



thin Ta films [2l|, and 
other InOx samples such as in reference [l7| ■ While many 
systems have shown scaling near the SIT with zv « 1.3 
[S 0, E3] ) no scaling consistent with what is expected 
for a true "dirty boson" system has been demonstrated. 
Figure 0] shows the scaling obtained for the highly dis- 
ordered sample of Fig. We note that this sample 
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FIG. 4: Top: Temperature scaling for sample "b" above. R c 
and H c are determined from the isotherms' crossing point. The 
resulting exponent zv is determined from the fit. The inset shows 
the crossing of the isotherms. Bottom: The same data as above 
but with scaling enforced with theoretical models for classical per- 
colation (zv = 4/3) and quantum percolation (zv = 7/3). 



exhibits very low H C /H C 2(0) and a critical resistance of 
R c w 6.5 kf2. Scaling near the crossing point was ob- 
tained in two steps. First, a best fit in which the expo- 
nent product zv was a free variable was employed. The 
scaling for each data set was based on its own actual 
crossing point. Following the procedure described previ- 
ously jjj, 1, 0, E3| ! the exponent product for this sample 
was found to be zv = 2.3±0.2 (see Fig.@Ji). We note that 
this product is what we normally get in scaling for the 
low resistance films, plus one, pointing towards a quan- 
tum percolation exponent. However, to test this result 
against specific theoretical predictions, we analyzed the 
data with two fixed values for zv predicted by percolation 
theory: zv — 4/3 for classical percolation and zv = 7/3 
for quantum percolation. It is evident from Fig. |4Jd that 
the higher exponent allows for significantly less spread of 
the data, reinforcing the free-fit result of 2.3 rather than 
a lower exponent. We have carried out a similar analysis 
on several other films of varying degrees of disorder 14 1 . 
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For other highly disordered films we found the scaling 
exponent to be consistent with v = 7/3. For the signif- 
icantly more weakly disordered samples (see or the 
sample in Fig. [5^) the scaling exponent was found to be 
v w 1.3 with a sometimes broadened crossing point and 
poorer scaling quality. Thus, the results in this regime 
essentially mimic those of MoGe. We have assembled all 
of our data, as well as various data from the literature, 
into the plot shown in Fig. [T] 

Let us finally return to a detailed discussion of Fig. [T] 
The vertical axis gives the conductance at the critical 
point in units of 4e 2 / h. The horizontal axis gives the field 
at the critical point, normalized by the estimated H C 2(0); 
the disorder decreases in the positive direction. For the 
InOx samples, H C 2(0) is calculated from the slope near 
T c q with some localization corrections as discussed in Q ; 
this field was also shown to be close to the peak of the 
magnetoresistance at low temperatures. For the MoGe 
samples H C 2 was measured from the low-field supercon- 
ducting transitions and the values are listed in references 
[1^ an d For the tantalum films, a discussion of the 
critical field and the coherence length is given in refer- 
ence 2l| . The error bars in Fig. [T] indicate the difficulty 
in determining H C 2 in the weakly insulating films. 

The data clearly display two branches. The samples 
in the top-right branch are very weak insulators with 
a critical conductance found to be well above the pre- 
dicted [6.45 kSl] -1 and a critical field in close proximity 
to the upper critical field H C 2- The scaling for all of 
these samples was found to be consistent with the clas- 
sical percolation model, with z = 1 and v w 1.35. The 
samples in the lower-left branch, on the other hand, are 
increasingly strong insulators with a critical conductance 
approximately the universal value, and scaling clearly 
shown in this paper to be consistent with v w 7/3. The 
samples near where the two branches meet (low G c but 
large H C /H C 2) show scaling that is more consistent with 
v = 4/3 than with 7/3 but of generally poor quality, 
which likely represents some crossover regime between 
the two branches. We therefore assert that there is a 
point (G*,H*) on our phase diagram which separates 
the tendencies towards insulating and metallic behavior, 
and we represent this point in Fig. [1] by x . 

In summary, as the disorder increases, the predictions 
of the dirty boson model and the expected quantum crit- 
ical behavior are clearly observed. Given the overall 
trend that the insulating strength increases as the dis- 
order increases, thereby suppressing any pair-breaking 
in the vicinity of the crossing point, it would appear 
that the assumption behind the dirty boson model — that 
the carriers are exclusively charge-2e bosons — is increas- 
ingly valid. We expect that in the limit of the most 
strongly disordered films, the predicted quantum critical 
behaviour should be observable down to lowest temper- 
ature. At the same time, the fact that the low disorder 
data clumps along a single branch maybe an indication 



of new physics in that regime. Such new physics 13j may 
be able to explain the widely observed metallic phase for 
high conductance films. 
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